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O^ I Abstract. Propp conjectured [15] that the number of lozenge tihngs of a semiregular hexagon 

f^ . of sides 2n — 1, 2n — 1 and 2n which contain the central unit rhombus is precisely one third of 

Cn ' the total number of lozenge tilings. Motivated by this, we consider the more general situation of 

^^ . a semiregular hexagon of sides a, a and b. We prove explicit formulas for the number of lozenge 

I tilings of these hexagons containing the central unit rhombus, and obtain Propp's conjecture 

Qs . as a corollary of our results. 

1. Introduction. Let a, b and c be positive integers, and consider a semiregular hexagon 
of sides a, b and c (i.e., all angles have 120 degrees and the sides have, in order, lengths a, 
/\ • bj c, a, 6, c). By a well-known bijection [2], the number of tilings of this hexagon by rhombi 

j^ ■ of unit edge-length and angles of 60 and 120 degrees (we call such a rhombus a lozenge and 

such tilings lozenge tilings) is equal to the number P(a, 6, c) of plane partitions contained in 
an a X 6 X c box. In turn, by a famous result of MacMahon [12], the latter is given by the 
product 

PC'M-flhfl'iiJTl^r (1-1) 
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2 M. CIUCU AND C. KRATTENTHALER 

The starting point of this paper is a conjecture of Propp [15] stating that for a semiregular 
hexagon of sides 2n — 1, 2n — 1 and 2n, precisely one third of its lozenge tilings contain the 
central lozenge. Call the lozenge tilings with this property centered. In this paper we 
consider the following more general problem: for a semiregular hexagon of sides a, a and 
6, how many of its tilings are centered? It is easy to see that such a hexagon has a central 
lozenge only if a and h have opposite parity. The two cases are addressed in Theorems 1 
and 2 below. 

For any nonnegative integer m and positive integer n define 



(2n)!2(2m)!(m + 2n-l)! /^ (-l)'^-*-i (m + n - 2)2^ 
^^"^^ ""' ~ 2-n!2m!(2m + 4n-2)! V ^ (2n - 2i - I) IX^ 

where the shifted factorial (a)^ is defined by {a)k '■= a{a + 1) ■ • ■ (a + /c — 1), k > 1, and 
(a)o := 1. 

Theorem 1. Let m be a nonnegative integer and n a positive integer. The number of 
centered lozenge tilings of a semiregular hexagon with sides 2n — 1, 2n — 1 and 2m is 
Q{m, n)P{2n - 1, 2n - 1, 2m). 

Theorem 2. Let m and n be positive integers. The number of centered lozenge tilings of a 
semiregular hexagon with sides 2n, 2n and 2m — 1 is Q{m^ n)P(2n, 2n, 2m — 1). 

In the case when m equals n, the expression Q{m,n) evaluates to 1/3 (this is due to a 
remarkable simplification of the sum in (1.2) in this case). Thus, the statement in Propp's 
conjecture follows from Theorem 1. 

Corollary 3. Let n be a positive integer. Exactly one third of the lozenge tilings of a 
semiregular hexagon with sides 2n — l, 2n — l and 2n contain the central lozenge. The same 
is true for a semiregular hexagon with sides 2n, 2n and 2n — 1. 

Theorems 1 and 2 have also been found (in an equivalent form) by Helfgott and Gessel 
[5, Theorem 2], using a completely different method. They also obtained Corollary 3 as a 
corollary of their results. 

For m ^ n the sum in (1.2) does not seem to simplify. However, if m and n approach 
infinity so that their ratio approaches some non- negative real number a, Q{m,n) turns out 
to approach the value - arcsin(l/(a + 1)). 

Corollary 4. Let a be any nonnegative real number. For m ~ an, the proportion of the 
rhombus tilings that contain the central rhombus in the total number of rhombus tilings of a 
semiregular hexagon with sides 2n — 1, 2n — 1 and 2m is ^ - arcsin(l/(a + 1)) as n tends 
to infinity. The same is true for a semiregular hexagon with sides 2n, 2n and 2m — 1. 

Remark. Using the bijection [2] between lozenge tilings and plane partitions, the statement 
of Theorem 1 can be interpreted as follows. Let V be the set of plane partitions (aij) of 
square shape (2n — 1)^"^"^, with entries between and 2m. Let Vk be the subset consisting 
of the plane partitions for which an+k,n+k = rn + /c, for — min(n — l,m) < k < min(n — 1, m). 
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Then the number of elements in the union of the VkS is Q{m, n)P{2n — 1, 2n — 1, 2?ti) (this 
union is clearly disjoint). 

A similar interpretation can be given to the statement of Theorem 2. 

The rest of the paper is devoted to giving proofs of Theorems 1 and 2, and of Corollaries 3 
and 4. In Section 2 we provide proofs of Corollaries 3 and 4, and outline the proofs of 
Theorems 1 and 2, the latter consisting of several steps. The details of these steps are 
then given in detail in the subsequent sections. These steps are the following. First, an 
application of the first author's Matchings Factorization Theorem [1, Theorem 1.2] allows to 
reduce our problem to the enumeration of lozenge tilings of simply-coimected regions. This is 
described in Section 3. Then, in Section 4, lozenge tilings are translated into nonintersecting 
lattice paths. By the main theorem on nonintersecting lattice paths [3, 4], the number (s) 
of nonintersecting lattice paths that we are interested in can be immediately written down 
in form of a determinant (see Lemmas 11, 12, 13). Finally, in Section 5, these determinants 
are evaluated (see Lemmas 15 and 16). 

2. Outline of proofs. Here we outline the proofs of Theorems 1 and 2 and we deduce 
Corollaries 3 and 4. We fill in the details in the subsequent sections. Denote by L{R) the 
number of lozenge tilings of the region R. 

Proof of Theorem 1. In Section 3 it is shown that the number of centered lozenge 
tilings of a semiregular hexagon with sides 2?i — 1, 2n — 1 and 2m equals 2^"~^ times the 
product of the number of lozenge tilings of two regions of the triangular lattice, if "*" and 
H~ (see (3.2) and Figure 3.3). Then, in Section 4 we use the Gessel-Viennot method of 
nonintersecting lattice paths to obtain determinantal expressions for L{H~^) (see Lemma 11) 
and L{H~) (see Lemma 13). Finally, in Section 5 we evaluate these determinants (see 
Lemmas 15 (with A^ = 2n — 2) and 16). After some manipulation of the expressions on the 
right-hand sides of (5.2) and (5.3) one obtains the statement of the Theorem. D 

Proof of Theorem 2. We proceed along the same lines as in the proof of Theorem 1. 
In Section 3 we show that the number of centered lozenge tilings of a semiregular hexagon 
with sides 2n, 2n and 2m — 1 equals 2^"~^ times the product of the number of lozenge tilings 
of two regions of the triangular lattice, H~^ and H~ (see (3.3) and Figure 3.4). The region 
H~ is the same as the one appearing in the proof of Theorem 1, while H^ differs from H~^ 
only in the sizes of its sides. The determinant evaluations in Lemmas 15 (with N = 2n and 
m replaced by m — 1) and 16 lead, after some manipulation of the expressions involved, to 
the statement of the Theorem. D 

Proof of Corollary 3. We have to compute the value of the expression on the right- 
hand side of (1.2) for m = n. Clearly, except for trivial manipulations, we will be done once 
we are able to evaluate the sum in (1.2) for m = n. 

We claim that 



E 



(_l)n-.-i (2n-z)2. _,„_inr=i'(6^-l)(6z + l) ^2.1) 



(2n-2z-l) iP (2?i-l)!!2 

(where the empty product is defined to be 1). Let us denote the sum by S{n) and its 
summand by F{n,i). We use the Gosper-Zeilberger algorithm [14, 19, 20] to obtain the 
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relation 

n{2n + lfF{n + 1, i) - 3n (6n - 1) (6?i + 1) F(n, i) = G{n, i + 1) - G{n, z), (2.2) 

with 

G{n, i) = (-1)''"' (-3 + 9 z - 6 z^ - 30 n + 62 m - 28 fn - 104 n^ + 104 m^ - 112 n^) 



X 



i^(2n-z + 2)2^-2 
(2?i-2z + l) i\^ 



Summation of the relation (2.2) from i = to i = n, little rearrangement, and division by n 
on both sides, leads to the recurrence 

(2n + lfS{n + 1) - 3 (6?i - 1) (6?i + 1) S{n) = 

for the sum in (2.1). (Paule and Schorn's [13] Mathematica implementation of the Gosper- 
Zeilberger algorithm, which is the one we used, gives this recurrence directly.) Since S{1) = 
1, and since the right-hand side of (2.1) satisfies the same recurrence, equation (2.1) is 
proved, and, thus, the Corollary also. D 

Proof of Corollary 4. First let a > 0. We have to determine the limit of Q{m, n) 
as n tends to infinity, and where the relation between m and n is fixed by to ~ an. Clearly, 
the "difficult" part of this asymptotic computation is to find the asymptotics of the sum 
in (1.2). It turns out that it is convenient to manipulate this sum first, before taking the 
limit n -^ oo. We reverse the order of summation in the sum in (1.2), and then are able to 
rewrite the sum using the standard hypergeometric notation 



r-*- s 



Oi, 



k=0 



{ai)k- ■ ■ {ar)k k 



k\{b 



1 k ■ ■ ■ 



{bs 



(2.3) 



(The reader should notice that the sum in (1.2) cannot be directly converted into a hyper- 
geometric series since the upper bound on the summation index, n — 1, is "artificial", i.e., 
the sum is changed if we extend the range of summation to all nonnegative i. If however 
the order of summation is reversed, this problem is removed because of the presence of the 
term il in the denominator of the summand.) Thus, expression (1.2) is converted into 



(2?i)!2 (2m)! {m + 2n-l)\{m + l)2n-2 
2-n!2m!(2m + 4n-2)! (n-1)!' 



-4-^3 



1, 2,1 -?^, 1 



n 



1 + m,2 — m — 2n 



3 ' 

' 2 



Next we apply Bailey's transformation for a balanced 4F3-series (see [17, (4.3.5.1)]), 



4-^3 



a, 6, c, —N 
e,/, l + a + 6 + c — e — / 



A^ 



;1 



(e-a)Af (/ -a )N 

(e)7v(/)7V 



"4-^3 



-N,a,l + a + c-e-f-N,l + a + b-e-f-N 



1+a+b+c 



f-N,l + a-e-N,l + a-f-N 
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where A^ is a nonnegative integer. Thus we obtain the expression 



{2n)l^ (2m)! (m + 2n - 1)! {m)n-i {m + n + V 



n-l 



2 ■ nl^ ml {2m + 4n - 2) 



{n-iy/ 



-4-^3 



L 2' 



1 — n, 1, 1, ^ +n 
2 — m — n, 1 + m + n 



Now we substitute m ^ an and perform the hmit n -^ oo. Using StirUng's formula, it is 
easy to determine the hmit for the quotient in front of the 4 Fa-series. It is 



2y/a(a + 2) 
7r(a + l)2 ■ 



(2.4) 



For the 4F3-series itself, we may exchange limit and summation, because of uniform conver- 
gence. This gives 



lim 4F3 



1 



n. 



1,1, ^+n 



L 2' 



2 — m — n, l + m-|-n 



;i 



— 2F1 



1,1 

3 

L 2 



1 



(a + 1) 



(2.5) 



(Recall that to ~ an.) Combining (2.4) and (2.5), and using the identity (see [16, p. 463, 
(133)]) 



2-ri 



1,1 

3 ;^ 
2 



arcsm V2 



V^(l-^) 



in (2.5), we obtain the desired limit -arcsin(l/(a + 1)). 

Finally we address the case a = 0, which means that to ~ 0. Then the proportion (1.2) 
is arbitrarily close to the expression which results from (1.2) when to ^ 0. In particular, 
according to Theorem 1, this expression gives the proportion of the rhombus tilings that 
contain the central rhombus in the total number of rhombus tilings of the semiregular 
hexagon with sides 2n — 1, 2n — 1 and 0. That proportion is simply 1 since, trivially, there 
is exactly one such rhombus tiling, and it does contain the central rhombus. The value of 1 
agrees with the claimed expression - arcsin(l/(a -(- 1)) evaluated at a = 0. D 

3. Reduction to simply-connected regions. One useful way to approach certain tiling 
enumeration problems is to biject them with nonintersecting lattice paths, and then use the 
Gessel-Viennot determinant theorem [3, 4]. This approach seems to be especially appropriate 
if the entries of the Gessel-Viennot matrix have a simple expression. In the case of the 
(2n — 1) X (2n — 1) x 2to semiregular hexagon with the central lozenge removed (whose 
tilings can clearly be identified with the centered tilings we are concerned with) this is not 
quite the case. However, one can get around this using the Factorization Theorem for perfect 
matchings presented in [1, Theorem 1.2]. 

Consider the tiling of the plane by unit equilateral triangles, illustrated in Figure 3.1. 
Define a region to be the union of finitely many such unit triangles. Suppose the region R 
is symmetric with respect to the horizontal symmetry axis /. Suppose further that the unit 
triangles of R crossed by / can be grouped in pairs such that the two triangles in a pair share 
an edge, forming a rhombic tile. Let Ti, . . . , T^ be these rhombi. Let P be the zig-zag lattice 
path that borders the tiles T^ on their upper boundary (see Figure 3.2). Define i?"*" and R~ 
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Figure 3.2 



H 

Figure 3.3 




Figure 3.4 



to be the pieces of R above and below P, respectively. Then the Factorization Theorem of 
[1] implies 



L{R) =2''L{R+)L*{R- 



(3.1) 



where L{R) is the number of lozenge tilings of R, and L*{R~) is the weighted count of the 
lozenge tilings of R~ assigning weight 2~* to a lozenge tiling containing i of the rhombi 
Ti, . . . , Tfc (a similar corollary of the Factorization Theorem is given in Remark 2.3 of [1] for 
the square lattice). 

Let H be the region obtained from a semiregular hexagon with side- lengths 2n — 1, 2n — 
1, 2to by removing the central lozenge (see Figure 3.3). By (3.1) we obtain 



L{H) = 22"-2L(if+)L*(if- 



(3.2) 



where the regions H^ and H are indicated in Figure 3.3 (H^ is obtained from the piece 
above the zig-zag path by removing 2m forced lozenges, which are indicated by a shading). 
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Figure 4.1 

Similarly, applying (3.1) to the hexagonal region H (illustrated in Figure 3.4) with side- 
lengths 2n, 2n, 2m — 1 and central lozenge removed we obtain 



L{H) 



-)2n-l 



L{H+)L*{H- 



(3.3) 



(Indeed, the region obtained from the bottom piece in Figure 3.4 after removing the forced 
tiles, which are again indicated by a shading, is the same as the region H~ in (3.2); this 
explains the last factor in (3.3); H~^ is shown in Figure 3.4). 

4. Enumeration of lozenge tilings and nonintersecting lattice paths. In this section 
we make use of the standard encoding of lozenge tilings in terms of non-intersecting lattice 
paths (Figure 4.1. illustrates this in the case of the region H^). Thus we transform the 
problem of enumerating lozenge tilings of the regions that arose in Section 2 into the problem 
of enumerating certain families of nonintersecting lattice paths. This allows us to derive 
determinantal formulas for the number of lozenge tilings we are interested in. 
Our next Lemma exhibits an expression for L(H~^). 

Lemma 11. We have 



L{H+) 



det 

l<i,j<2n-2 



2n + m — i - 
m + i — j 



1 



(4.1) 



Proof. Using the correspondence with nonintersecting lattice paths, one obtains that the 
number of lozenge tilings of the upper "half hexagon" H~^ (see Figure 3.3) equals the num- 
ber of families (Pi, P2, ■ ■ ■ , P2n-2) of nonintersecting lattice paths (consisting of horizontal 
positive unit steps and vertical negative unit steps) in which Pi runs from Ai = {2i,i + m) 
to Ei = {2n + z — 1, i), i = 1, 2, . . . , 2n — 2. By the main theorem of nonintersecting lattice 
paths [4, Cor. 2; 18, Theorem 1.2], this number is given by the determinant 

det {\V{A,^E,)\), 

l<.t,j<2n — 2 

where V{A -^ E) denotes the set of paths from A io E and \M.\ is the cardinality of the 
set M.. Obviously, the number \V{Ai -^ Ej)\ of paths from Ai to Ej equals the binomial 



) . This establishes the Lemma. D 



2n+m — i- 
V m-\-i—j 

In the same way, we may derive a determinantal formula for L{H^). As we have already 
noted in Section 2, the only difference between the regions H~^ and H^ is in their side 
lengths. To be precise, the vertical sides of H^ have length m, while those of H^ have 
length TO — 1. On the other hand, the slanted sides of H~^ have length 2n — 2, while those of 
H^ have length 2n. Hence, we will get a formula for L{H~^) by replacing m by to — 1 and 
n by n + 1 in the formula for L{H~^). 
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Lemma 12. We have 

L(H^)= det ff'" + "^r\')V (4.2) 

D 

Next we turn to the number L{H~). In analogy to the preceding, we derive also a deter- 
minantal expression for L{H~), using nonintersecting lattice paths. However, the resulting 
determinant is not as "regular" as the preceding one, and is therefore harder to evaluate. 

Lemma 13. We have 

L(H-)= det ( (2n + m-z-l)! ((n + m - ^) tfty^n)\ 

i<i,j<2n-i\{m + i-jy.{2n-2i + jy.\ j if i = n j J ^ ' 

Proof. Using the correspondence with nonintersecting lattice paths, the number of lozenge 
tilings of the lower "half hexagon" H~ (see Figure 3.3) equals the number of families 
(Pi, P27 • • • 7 -P2n-i) of nonintersecting lattice paths (consisting of horizontal positive unit 
steps and vertical negative unit steps) in which Pi runs from Ai = {2i,i + m) to Ei = 
(2n + i, i), i = 1, 2, . . . , n — 1, n + 1, . . . , 2n — 1, whereas Pn runs from A^ — (2n + 1, n + to) 
to En — (3n, n) (i.e., it is the starting point A^ which deviates slightly from the "general" 
rule). In this count, horizontal steps originating from any Ai, i ^ n, count with weight 1/2. 
Again, by the main theorem of nonintersecting lattice paths, this number is given by the 
determinant 

det (\V{A^E,)Q, 



where |.|^„ denotes this weighted count. It is not difficult to see that the weighted count 

i+m— j/2) (2n+m,— i — 1 
{m+i-j)\ (2n-2i+j)! 



\V{Ai -^ Ej)\w of paths from Ai to Ej equals '-""t^, /-,^?! (2n-2z+'7)!^'" if ^ 7^ ^ and 



rn+m \ _ J (n "^ )■ otherwise. This establishes the Lemma. D 

Vm+n— J/ (m+n—j)'.j'. 

5. Determinant evaluations. Here we evaluate the determinants which appear in Lem- 
mas 11, 12 and 13. 

In order to compute the determinants in (4.1) and (4.2), we utilize the following deter- 
minant lemma from [6, Lemma 2.2]. 

Lemma 14. Let Xi, . . . ^Xn, A2j . . . , An, and B2j . . .Bn be indeterminates. Then there 
holds 



det Ux, + An)--- {X, + A,+i)(X, + B,) ■ ■ ■ (X, + B2] 

l<i,i<N \ 



n (x^-xj) n (^^-^o- (5-1) 

l<i<j<N 2<i<j<Ar 



D 



Now we can state, and prove, the evaluations of the determinants in (4.1) and (4.2), in a 
unified fashion. 
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Lemma 15. For any positive integer n there holds 



N 



N + m-i + l\\ y-f {N + m-i + l)\{i- l)!(2m + z + l)i 



i<?i<ivVV m + i-j ) ) ~ ^^ {m + %-l)\{2N-2i + l)\ ' ^ ' 



i=l 



Proof. We take (A^ + m - z + l)!/((m + z - 1)! (2A^ - 2z + 1)!) out of the i-th row of the 
determinant in (5.2), z = 1, 2, . . . , A^. Thus we obtain 

^ (iV + m-i + 1)! 



n 



^J^ {m + i- 1)! (2A^ - 2z + 1)! 

X det ((m + z- j + l)---(m + i-l)(A^-2z+j + 2)---(2A^-2z + l)). 

l<i,j<N^ ^ / \ V 

Taking (— 2)^"-^ out of the j'-th column, j = 1, 2, . . . , A^, we may write this as 



iv^JX- {N + m-i + 1)1 



(-2)(;)n 



fj^ {m + i- 1)! (2A^ - 2z + 1)! 

Now Lemma 14 can be apphed with Xi = i, Aj = —{N + j + l)/2, Bj = m — j + 1. After 
some simphfication one arrives at the right-hand side of (5.2). D 

The determinant in (4.3) evaluates as follows. 

Lemma 16. For any positive integer n there holds 

(2n + m-z-l)! j(^n + m-^) if i ^ n 



(let ^ ■ - < ^ 2 

i<i,j<2n-i\{m + i — j)\{2n — 2i+ j)\\ j if i = n 

n 2n-l /o I 'I'll 2n-2 



i=0 



(2?i-2z-l) i!^ 



Proof. The method that we use for this proof is also applied successfully in [10, 7, 8, 9, 11] 
(see in particular the tutorial description in [9, Sec. 2]). 

First of all, as in the proof of Lemma 15, we take appropriate factors out of the determi- 
nant. To be precise, we take (2n + m — i — l)!/((rM + i — 1)1 (4n — 2i — 1)!) out of the i-th 
row of the determinant in (5.3), z = 1, 2, . . . , 2n — 1. Thus we obtain 

T-r (2n + m — i — 1)1 



J- (m + z-l)!(4 



n-2z-l)! 

X det ( (m + z- 7' + l)---(m-F-z- l)(2?i-2z + ?' + l)---(4n-2z- 11 

(n + m — |) if z 7^ n 
J if z = n 

(5.4) 
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Let us denote the determinant in (5.4) by D{m; n). Using the notation of shifted factorials, 
this means that 

D{m;n):= det ( {m+i-j + l),-i {2n-2i+j + l)2n-j-i I ^""^ "^ " ^^ '^'^'']]. 

(5.5) 
Comparison of (5.3) and (5.4) yields that (5.3) will be proved once we are able to establish 
the determinant evaluation 

TT'^ ('97 - 1^2 "~^ 
D{m; n) = 22(-i)('^-2) ^^^='^_'^- H ((^ + ^)2-2. {m + i + l/2)_, ) 

[n i). 



i=l 

g (-ir-Mm + n-.),. 
^(2n-2i-l) i!2 ^ ' 

i=0 ^ ' 

For the proof of (5.6) we proceed in several steps. An outline is as follows. In the first 
step we show that HiLi (^ + 02n-2i is a factor of D{m] n) as a polynomial in m. In the 
second step we show that HlLi {'m' + i + 1/2) n-i is a factor of 0(771; n). In the third step 
we determine the maximal degree of D{m; n) as a polynomial in m, which turns out to be 
(2n + l)(n — 1). From a combination of these three steps we are forced to conclude that 

n-l 

D{m; n) = Yl {{m + i)2u-2i {m + i + l/2)^_i )P(m; n), (5.7) 

where P{m; n) is a polynomial in m of degree at most 2n — 2. Then, in the fourth step 
we show that P{m;n) = P{1 — 2n — m;n). And, in the fifth step, we evaluate P{m;n) at 
m = 0, —1, . . . , — n + 1. Namely, for to = 0, —1, . . . , — n + 1 we show that 

P{m;n) = (_i)" 2^("-^)(-^) ^L^ ^^\~/^'' ^^ ^^^^'^ (5.8) 

Clearly the latter two properties determine a polynomial of maximal degree 2n — 2 uniquely. 
As is easy to check, the sum in (5.6) has the first property, too, namely that it is invariant 
under replacement of m by 1 — 2n — m. Since in the sixth step we prove that for m = 
0, —1, . . . , —n + 1 we also have 

n-l ^_^^^_,_i ^^^^_--^^^ ^ (^+n)m+n-l 

^(2n-2i-l) z!2 ^ > 2{^-nU+^ 

we are forced to conclude that 

P(^; n) = 2^(-^)(-'> n-^i (2. - 1)P g (-1)"-' (m + n- .).. 

^ ' ^ (n-l)! ^(2n-2z-l) iP ^ ' 
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This would finish the proof of the Lemma since a combination of (5.7) and (5.9) gives (5.6), 
and thus (5.3), as we aheady noted. 

Step 1. nr=i (^ + 02n-2i is a factor of D{m; n). For i between 1 and n — 1 let us consider 
row 2n — i of the determinant D{m;n). Recalling the definition (5.5) of D{m;n), we see 
that the j-th entry in this row has the form 

{m + 2n- i - j + l)j-i (-2n + 2i + j + l)2n-j-i in + m-- 

Since (-2n + 2i + j + l)2n-j-i = for j = 1, 2, . . . , 2n - 2z - 1, the first 2n - 2z - 1 
entries in this row vanish. Therefore (m + i)2n-2i is a factor of each entry in row 2n — i, 
i = 1, 2, . . . , 71 — 1. Hence, the complete product YYi=i ("^ + ^)2n-2i divides 0(771; n). 

Step 2. nr=i (^ + ^ + l/2)n-i is a factor of D{m;n). Let us concentrate on a typical 
factor (m + j + I + 1/2), l<j<n — l,0<l<n — 2. We claim that for each such factor 
there is a linear combination of the rows that vanishes if the factor vanishes. More precisely, 
we claim that for any j, I with l<j<n — l,0</<n — 2 there holds 



(-1) 



j_i (n-j-l- l)j 



4j^n-j-l), 

■ (column (2n - 2j -21-1) of D{-j - I - 1/2; n)) 

2n-2l-l , . s 

+ E ( ^-^9/ 9 ^J-(column.ofD(-j-/-l/2;n)) = (5.10) 

^ V'5 + J + 2f - 2n + 1/ 



s=2n-j-2Z 

if j + / < n, and 



2i+i 



L^ \ s + l-n ) ■ (col^"^^ '5 of ^(-^ ~^~ 1/2;^)) 



s^n- 



(-l)""^(j + /-n + 1)^-^-1 

+ ^n-l-l ^j + I _ n + 2)n-l-l 
2j+2l-2n+3 
X 

s=l 



G 1 ^ ' 



if J + / > n. 

In order to verify (5.10), for j + / < n we have to check 

2n-2Z-l / 

.^2n^-2.-A^+^' + 2/-2n + l 



■ (n - J - / - I - -) (z - J - / - s + i),_i (2n - 2i + s + 1; 



2n-s-l 



0, 
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which is (5.10) restricted to the i-th row, i ^ n (note that the entry in column 2n — 2j — 2/ — 1 
of D{—j — I — 1/2; n) vanishes in such a row), and 

j_l {j + l-n+ |)2n-j-2Z-2 (2n - 2j -21- l)2j+2Z+l 

2n-2l-l , . X 

+ y] ( , • , o"? O , . )(^-i-^-S+i)s-l(s)2n-s = 0, 

\s + ? + 2/ — 2?i + 1/ ^ 

which is (5.10) restricted to the n-th row. Equivalently, using the standard hypergeometric 
notation 

ai, . . . ,ar 






6i,...,6s 



A;=0 



(ai)fc ■ ■ ■ {ar)k k 
kl {bi)k---{bs)k 



this means to check 



{i + l-2n+ l)2n-j-2l- 



.2 (4n -2i-j- 2l)j+2i 3F2 



1+ j, -i -i-l + 2n,-j 
j, -2i -j-2l + 4n ' 



(5.12) 



and 



j_l (j + / - n + |)2n-j-2Z-2 (2?1 - 2j - 2/ - l)2j+2Z+l 

^ ^ 4^(n-J-/), 



+ (/ - n + |)2n-J-2Z-2 (2n -j -21- l)j+2Z+l 2i^l 



/ + n, -j 



2 - ' ■"7 J . 1 

-l-j-2/ + 2n' 



0. (5.13) 



In order to verify (5.11), for j + / > n we have to check 



-I 



s + I — n 



2 2 



{i - j - I - s + l)s-i {I - 2i + 2n + s] 



2n-s-l 



+ 



(-l)""'(i + /-n + f)„-z-i 



2j+2l-2n+3 



s=l 



2j + 2l-2n + 2 
s-1 



n-j-l- 



s 1 



2 ^]{i-j-l-s+l)s-iil-2i + 2n + s] 
which is (5.11) restricted to the i-th row, i ^ n, and 



2n— s — 1; 



2j+l 



2j + / - 71 + 1 



{n- j -I- s + l)s-i (s)i 



S + I — n ' ^ 



\n — l 



+ 



(_l)--(j + /_n + f)._z_ 



4«-'-i (j+/-n + 2)^_;_i 

2j+2Z-2n+3 



s = l 



2j + 2/ - 2n + 2 



s-1 



{n-j -I- s + \)s-i{s) 



2n—s: 
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which is (5.11) restricted to the n-th row. Equivalently, using hypergeometric notation, this 
means to check 



(n — 2j — / — 1) 1 

7^ (2 +^- J - n)n-i-i (3n - 22 - / + l)n+l-l 



X 3^2 



2 + 2j + / - 71 i - z + j + n, -1 - 2j - / + n 



1 + 2j + / - n, 1 - 2z - / + 3?i 



;i 



+ 



(_l)--'(n-j-/-l)(|+j+/-n) 



-i-i 



{2n - 2i + 2)2n-2 



X 3F2 
and 

(i - j)n-l-l {n - l)n+l 2F1 



4"-^-^ {2 + j + I - n)n-i-i 

3 + 2j + 2/ - 2n, I - z + j + /, -2 - 2j - 21 + 2n 



2 + 2j + 21- 2n, 2-2i + 2n 



;1 



(5.14) 



l+j,-l-2j-l + n 
-l + n 
-I /1^ /3 



+ 



(-l)""'(l)2n-l(f +i + /-?^)n-Z-l 



-l-l 



X2F1 



{2 + j + I ~ n)^-i-i 
I + J + / - n, -2 - 2j - 2/ + 2n _ 



1 



;i 



0. (5.15) 



We start with the proof of (5.12). We apply the contiguous relation 



3-^2 



a, ^1,^2 
-61,-62 



3-^2 



a — 1, Ai, A-2 
-61,-62 






a,Ai + 1,^2 + 1 

Si + 1,-82 + 1 



to the 3F2-series in (5.12). Since in this case a — 1 = Bi, parameters cancel inside the two 
3F2-series on the right hand side of the contiguous relation, leaving two 2-^1-series instead. 
Thus, (5.12) is turned into 



{i + l-2n+^)2n-j-2l- 



.2{4n-2i-j-2l)j+2ii2Fi 
(2n-z-/-i) 



4n - 2i - j - 2/ 



i^i 



2n — i — I — |, — i . -, 
4n - 2z - j - 2/ ' 

^ — i — I + 2n,l— j 
1 - 2i - j - 21 + 4n' 



0. 



Each of the two 2-P'i-series can be evaluated by means of the Chu-Vandermonde summation 

(see [17, (1.7.7); Appendix (III.4)]), 



iFi 



a,-N 



1 



(c — a) 



N 



(c) 



N 



(5.16) 



where N is a, nonnegative integer, and thus (5.12) follows upon minor simplification (the 
terms in big parentheses cancel each other). 

To the 2-P'i-series in (5.13) Chu-Vandermonde summation (5.16) can be applied directly, 
and it yields the desired result. 
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The verifications of (5.14) and (5.15) are similar. The reader wiU have no difficulties to 
fill in the details. 

This finishes the proof that the product nr=i ("^ + ^ + l/2)n-i divides D{m;n). 

Step 3. D{m]n) is a polynomial in m of maximal degree (2n + l)(?i — 1). Obviously, 
the degree in m of the (z,j)-entry in the determinant D{m;n) is j for i ^ n, while it is 
j — 1 for i = n. Hence, in the defining expansion of the determinant, each term has degree 

(E-:7'j)-i = ft^)-i = (2^+i)(^-i)- 

Step 4- P{rn; n) = P(l — 2n — m; n) . We claim that there holds the relation 



D{m; n) = (-1)""'D(1 - 2n - m; n). 



(5.17) 



It is clear from the definition (5.7) of P{ni;n) that (5.17) immediately implies the desired 
relation P(m; n) = P(l — 2n — m;n). 

We prove (5.17) by, up to sign, transforming the determinant D{ni;n) into the determi- 
nant -D(l — 2?i — to; n) by a sequence of elementary column operations (which, of course, 
leave the value of the determinant invariant). To be precise, for j = 2n — 1, 2n — 2, . . . , 2, 
we add 



E 

fc=i 



J -I 
k-1 



(column k of D{m; n)) 



to column j. Thus, in the new determinant, Di{'m;n) say, the (i, j)-entry is 



E 



fc=i 
for i ^ n, and 



\[m + n - -\{m -{- i - k + l)k-i {2n - 2i + k + l)2n-k-i 



E 

k=i 



J-1 
k-1 



{m + n-k + l)fc_i {k)2n-k 



for i = n. Using hypergeometric notation, the (i, j)-entry of Di{m; n) is 



m + n--]{2-2i + 2n)2n-2 3F2 



2 — 2m — 2n, 1 — i — ni,l — j 
1 - 2m - 2n, 2 - 2i + 2n 



(5.18) 



for i j^ n, and 



(l)2n-i 2F1 



1 — m — n, 1 — J 



(5.19) 



for i = n. The 3F2-series in (5.18) can be evaluated in the same way as we evaluated 
the 3F2-series in (5.12). The 2-f^i-series in (5.19) is easily evaluated by Chu-Vandermonde 
summation (5.16). Thus, we obtain that the (i,j)-entry of Di{ni;n) is given by 



(_1)J l^i-±^jn-n){2 + i-j-m- 2n)j-i {1 - 2i + j + 2n)2n-j-i (5.20) 



THE NUMBER OF CENTERED LOZENGE TILINGS 15 

for i ^ n^ and by 

(-l)^-'i (2 - J - m - n),_i (j + l)2n-,-i (5.21) 

for i = n. Now, the expressfon (5.20) is exactly (—1)-' times the (i, j)-entry of D{1 — 2n — 
m;n), while the expression (5.21) is exactly (— 1)-^"^ times the (n, j)-entry of D(l— 2n— ?B;n). 
Hence, relation (5.17) follows immediately, implying P{m;n) = P{1 — 2n — m;n), as we 
already noted. 

Step 5. Evaluation of P{m; n) at m = 0, —1, . . . , —n + 1. The polynomial P(m; n) is 
defined by means of (5.7), 

n-l 

D{m; n) = n ((^ + ^)2n-2^ (m + z + l/2)^_i )P(m; n). (5.22) 

So, what we would like to do is to set m = — e, e being one of 0, 1, . . . , n — 1, evaluate 
D{~e;n), divide both sides of (5.22) by the product on the right-hand side of (5.22), and 
get the evaluation of P{m; n) aX m = —e. However, the product on the right-hand side of 
(5.22) unfortunately (usually) is zero for m = — e, < e < n — 1. Therefore we have to find 
a way around this difficulty. 

Fix an e with < e < n — 1. Before setting m = — e in (5.22), we have to cancel (m + eY 
on the right-hand side of (5.22). To accomplish this, we have to "generate" these factors on 
the left-hand side. We do this by adding 

2e+2j-2n 



y \ 1 ■ ] ■ (column /c of D (to; n)) 



k=j+i 

to column j, j = 2n — 2e, 2n — 2e + 1, . . . , 2n — e — 1. Thus, in the new determinant the 
entry in the i-th row in such a column is 



(m + n- -jl _. Um + t-k + l)fc_i {2n - 2t + k + 1)2 

k=j 

a i j^ n, and 

Ef2e + j-2n\. , , ,,, 

I , . j{m + n-k + l)k-i{k)2n-k 



k=j 

ii i — n. In hypergeometric terms this is 



(m + n - -]{m + i - j + l)j_i {2n - 2i + j + l)2n-j-i 

1 + J — 2m — 2n, —i+j — m, — 2e — j + 2n 



X 3F2 
if i 7^ n, and 

{m + n- j + l)j-i {j)2n-j 2F1 



■ 1 
j -2m- 2n, l-2i+j + 2n ' 



j — m — n, —2e — j + 2n 
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a i = n. Again, the 3F2-series can be evaluated in the same way as before the 3-F2-series 
in (5.12), while the 2-^1-series is evaluated by means of the Chu-Vandermonde summation 
(5.16). Thus, we obtain that the (i,j)-entry, 2n — 2e < j < 2n — e — 1, of the modified 
determinant is given by 

(m + e) {m + i - j + l)j-i (m + 2n - i)2e+j-2n (2e - 2i + 2j + l)^ri-2e-2j-i 

if i 7^ n, and by 

(m + e) {m + n-j + l)e+j-n-i (m + e + l)e+j-n-i (2e + 2j - 2n)4„_2e-2j 

ii i = n. Clearly, {m + e) is a factor of each entry in the j-th column of the modified 
determinant, 2n — 2e < j < 2?i — e — 1. Therefore, we may take (m + e) out of the j'-th 
column, j = 2n — 2e, 2n — 2e + 1, . . . , 2n — e — 1. The remaining determinant, D2{m; n) say, 
is then defined as D2{m;n) = deti<ij<2n-i(-E'ij), where for 1 < j < 2n — 2e — 1 and for 
2n — e < J < 2?i — 1 the entry Eij is given by 

i?., = (m + z-j + l),_i(2n-2z+j + l)2.-,-i|(^ + 'J-2) ![ ^ ^ |^ (5.23) 

(i.e., Eij equals the (i, j)-entry of D(m; n), as given by (5.5), in that case), and for 2n — 2e < 
j < 2n — e — 1 the entry Eij is given by 

({m + i-j + l)j_i (m + 2n - i)2e+j-2n (2e - 2z + 2j + l)^^_2e-2i-i i^i^n 
Eij = \ 

[ {m + n~j + l)e+j-n-i {e + m + l)e+j-n-i (2e + 2j - 2n)4n-2e-2j ifi = n. 

(5.24) 
Due to the manipulations that we did, the new determinant D2{m]n) is related to the 
original determinant D{m] n) by 

D{m; n) = {m + e)^ D2{'m; n). 

Substituting this into (5.22), and rearranging terms, we get 

e 

P(m; n) = D2{m; n) J]^ ((m + i)e-i (m + e + l)2n-i-e-i {m + i + \I2)^_^ ) 

n-l 

J] ((m + i)2n-2^(m + z + l/2)^_J"\ (5.25) 



n-l 
X 

i=e+l 



Now we may safely set ra = —e. So, what we need in order to obtain the evaluation of 
P{m; n) 3it m = — e is the evaluation of the determinant D2{—e] n). 

In order to determine the evaluation of 1^2 (~e; n), we observe that D2{—e; n) has a block 
form which is sketched in Figure 10. 
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J = l 



3 = 
2n — e 



* 


* 


* 


* \ 


* 








X 






i^ 1 

i = e + l 
z = 2?i — e 



2n 



2n - 2e 



3 



2n 



The block form of D2{—e;n) 
Figure 10 

The figure has to be read according to the following convention: If a block is bounded by 
horizontal lines marked as i = hi and i = /i2 and vertical lines marked as j = vi and j = V2, 
then the block consists of the entries that are in rows i = hi, hi + 1, . . . , h2 — 1 and columns 
j = vi,vi + 1, . . . ,V2 ~ I- It is an easy task to check from the definitions (5.23) and (5.24) 
of the entries of D2{m;n) that indeed the lower left block of D2{—e;n), consisting of the 
entries in rows z = 2n — e, 2n — e + 1, . . . , 2n — 1 and columns j = 1, 2, . . . , 2?i — 2e — 1, and 
the lower right block, consisting of the entries in rows i = e + l,e + 2, ...,2?i — 1 and columns 
j = 2n — e, 2n — e + 1, . . . , 2n — 1 are blocks of zeroes. Hence, the determinant D2{—e; n) 
factors into the product 

det(Si) det(S2) det(S3), 

where Bi is the middle left block, consisting of the entries in rows i = e+1, e+2, . . . , 2?i — e — 1 
and columns j = 1, 2, . . . , 2n — 2e — 1, where B2 is the lower middle block, consisting of the 
entries in rows i = 2n — e, 2n — e + 1, . . . , 2?i — 1 and columns j = 2n — 2e, 2n — 2e + l, . . . , 2n — 
e — 1, and where S3 is the upper right block, consisting of the entries in rows i = 1,2, ... ,e 



and columns j = 2n 



e, 2n 



e + 1. 



2n-l. 



As is indicated in Figure 10, the blocks Bi and B2 are lower and upper triangular matrices, 
respectively. Hence, their determinants are easily computed. 
The determinant of the block S3 is the determinant 



det {{-2n + i- j + 2)2n-e+j-2 {4n - e - 2i + j] 

l<*J<e 



1 



J 



)• 



(All the entries of S3 are given by the first case of formula (5.23) because e < n - 
We take {—2n + i + l)2n-e-i out of the i-th row of this determinant, i = 1, 2, . . . , e, 
(— 2)*^"-^ ( ~2~'' ) out of the j-th column, j = 1, 2, . . . , e. Thus we obtain 



and 



(-2)«)n 



1 — e — i 



{-2n + i + l)2n-e-l 



i=l 



X det ( (i — 2n + - 

l<i,j<e \\ 2 



i-2n + 



e-J 



(i - 2n - i + 2) ■ ■ ■ (z - 2n) 
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The latter determinant is easily evaluated using Lemma 14 with N — e, Xi = i, Aj — 
-2n + {e-j + l)/2, Bj = -2n -j + 2. 

This finishes the desired evaluation of D2{—e;n), and, via (5.25), of P{—e;n) for e = 
0, 1, . . . , n — 1. If everything is put together and simplified, the result is exactly (5.8). 

Step 6. Evaluation of the sum on the right-hand side of (5.3) at m = 0, —1, . . . , —n + 1. 
We claim that ifO<e<n — Iwe have 



"' (-1)'^-'-^ in-e-^)2^_, ,,ni^ + n 



^^^ (2n - 2z - 1) ^^ - ^ '^ 2(i-n)._e- ^^''^^ 



This is seen by first rewriting the sum in (5.26) as 



(_1)- - i(i_^)^(i + e-n)an-e), 



Since < e < n — 1, the term (1 + e — n)i will make the sum terminate at i = n — e — 1, and 
so at i = n — 1 latest. Therefore we may extend the range of summation to all nonnegative 
numbers z, without altering the sum. Then we can write the sum in hypergeo metric notation 
as 



"3-^2 



l-2n 



^— n, n — e,l + e — n 
lA-n ' 



"' 2 

This 3F2-series can be evaluated by means of the Pfaff-Saalschiitz summation (see [17, 
(2.3.1.3); Appendix (III.2)]), 



3-^2 



a, 6, — A^ _ ^ 

c,l + a + b-c- N' 



(c-a)jv(c-6)Ar 
(c)jv(c-a-6)Ar' 



where A^ is a nonnegative integer. Thus we arrive at the right-hand side of (5.26). 
This completes the proof of the Lemma. D 
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